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Abstract 


Within a permanently shadowed lunar crater the horizontal flow of solar 
wind is obstructed by upstream topography, forming a plasma wake that 
electrostatically diverts ions toward the crater floor and generates a surface 
potential that can reach kilovolts. In the present work kinetic plasma simulations 
are employed to investigate the morphology of a lunar crater wake during passage 
of a solar storm. Results are cast in terms of leading dimensionless ratios including 
the ion Mach number, ratio of crater depth to plasma Debye length, peak secondary 
electron yield, and electron temperature vs. electron impact energy at peak 
secondary yield. This small set of ratios allows generalization to a much wider 
range of scenarios. The kinetic simulation results are fed forward into an 
equivalent- circuit model of a roving astronaut. In very low-plasma-current 
environments triboelectric charging of the astronaut suit becomes effectively 
perpetual, representing a critical engineering concern for roving within shadowed 
lunar regions. Finally, simulated ion fluxes are used to explore sputtering and 
implantation processes within an idealized crater. It is suggested that the physics of 
plasma mini-wakes formed in the vicinity of permanently shadowed topography 
may play a critical role in modulating the enigmatic spatial distribution of volatiles 
at the lunar poles. 

1. Introduction 

From a human exploration perspective permanently shadowed lunar craters 
have the advantage of being thermally stable, potentially volatile-rich, [Paige et al. 
2010] and in some locations providing nearby accessibility to steady sunlight as an 



energy source. However, within a shadowed crater the direct flow of solar wind is 
obstructed by the upstream crater wall, forming a complicated plasma wake that 
modulates the particle fluxes and electrostatic environment far downstream 
[Zimmerman et al. 2011; Farrell et al. 2010]. Astronauts and robots venturing onto 
the crater floor will accumulate charge through direct exposure to spatially 
inhomogeneous plasma and secondary currents, and will undergo triboelectric 
charging due to frictional contact with the surface [Jackson et al. 2011], These 
charging processes potentially represent an electrical discharge hazard to the 
astronaut suit as well as other infrastructure including sensitive instrumentation 
[Jackson et al. 2011]. Furthermore, the processes of ion implantation and sputtering 
are thought to play a role in the local balance of volatiles, which may someday be 
exploited to produce water and other useful compounds in situ [Shearer and Neal, 
2011; Neal, 2009], Understanding the environments within permanently shadowed 
regions (PSRs) of the Moon, under the wide variety of solar wind conditions 
experienced at 1 AU, will feed forward into future in situ mission planning and 
provides a unique laboratory for exploring the physics of plasma wakes. 

The overarching goal of this work is to characterize the crater wake plasma 
environment, its effects on astronaut charging, and the volatile balance within an 
idealized polar crater during passage of a coronal mass ejection [CME]. This is 
achieved by coupling a kinetic plasma model of wake structure [Zimmerman et al. 
2011] with a numerical model of object charging [Jackson et al. 2010] and 
experimental and theoretical results on volatile sputtering [Johnson et al. 1989; cf. 
Killen et al. 2012 and references therein]. It will be shown that the coupled 



wake/surface system is highly sensitive to fluctuations in solar wind conditions 
during the various stages of the passing CME, and the charge state of a roving 
astronaut will vary dramatically within the respective conditions of each stage. 
Furthermore, while previous work suggested that wake physics within polar craters 
diverts solar wind protons toward the crater floor [where they may represent a 
significant source of implanted hydrogen) [Zimmerman et al. 2011, Farrell et al. 
2010], one stage of the CME introduces a significant population of heavier, highly- 
charged He ++ ions that may also perturb the balance of volatiles at the surface. 

For convenience a CME event observed by Skoug et al. 1999 is divided into 
four distinct temporal stages: typical solar wind; a dense, hot shock; a faster- 
flowing, cooler early cloud; and, a dense, cold late cloud (e.g., Table 1) [Farrell et al., 
submitted to this issue, 2012]. The constituent environmental conditions are 
simulated using the open-source kinetic plasma code XOOP1C [Verboncoeur et al. 
1995] incorporating an idealized, dielectric lunar surface and open boundary 
conditions on the solar wind plasma. Details of the kinetic code and simulation set, 
as well as conclusions regarding variations in wake structure during passage of the 
storm are presented in Section 2. Particle currents are extracted at representative 
locations at the surface and fed forward into the astronaut charging model of 
Jackson et al. 2011, modified to include an accounting of secondary electron 
emission from the astronaut suit. Section 3 contains a description of this charging 
model and presents results on the charging characteristics of insulating and 
conducting astronaut suits during the storm. Finally, Section 4 contains a discussion 



of the implications for volatile sputtering and sequestration due to the various 
proton and heavy ion fluxes encountered during the storm. 

2. Kinetic simulations of wake density, electric field, potential, and 
particle fluxes 

2.1. Kinetic simulation setup 

The simulation approach is largely as described in Zimmerman et al., 2011. 
The open-source 2D kinetic plasma code XOOPIC [Verboncoeur 1995] is used to 
simulate an idealized, step-like, dielectric crater topography. This simplification of 
shape, which is illustrated in Fig. 1, reveals the governing physics in a clear and 
straightforward manner. Maxwellian populations of electrons and protons are 
injected at the upstream and overhead boundaries to provide a continuous flow of 
plasma particles into the system. Slight modifications were made to the code to 
enable a self-consistent flow of particles through the downstream boundary. 
Electrons with horizontal velocity greater than 2 Vsw, where v*w is the solar wind drift 
speed (measured in m s _1 ) are recycled at the downstream boundary to simulate the 
backward-moving portion of the electron population as it responds to upstream 
conditions. All other particles exiting through open boundaries of the simulation 
domain are discarded. Otherwise, the particle trajectories evolve according to the 
locally self-consistent electric field including contributions from charges that have 
impacted and been recorded at the surface. 

The ion population is represented using a downward-moving, half- 
Maxwellian velocity distribution to capture the fact that upward-moving ions will 
have been absorbed by the dayside lunar surface far upstream. The magnetic field is 



assumed to be zero everywhere, as a typical electron undergoes about 3% of a 
gyration (about an interplanetary magnetic field line) while traversing the entire 
depth of the crater; the validity of this assumption has been verified through a 
handful of simulations with a nonzero interplanetary magnetic field. It is known 
that surface-reflected protons from the lunar dayside can reach polar and nightside 
regions after gyration within the interplanetary magnetic field [Nishino et al, 2009]; 
however, non-ambipolar ion entry into polar craters is beyond the scope of the 
present paper, and we do not attempt to model global-scale ion pickup in the 
present, regional-scale simulations. 

The electric field components at the crater surface are affected by the locally 
accumulated surface charge according to the boundary conditions imposed by 
Maxwell's equations [Verboncoeur et al. 1995], and the surface charge is naturally 
zero at the open solar wind boundaries. An electric potential gauge of zero volts is 
imposed at the overhead boundary. Secondary electrons are emitted by primary 
electron impacts using the yield-vs.-impact energy curve of Vaughan 1989, 
parameterized to model the experimentally-measured yield of lunar regolith [Willis 
1973]. Secondary emission by ion impacts (i.e. sputtering) is not included in the 
present model. Small variations in the net ion current due to sputtering of other ions 
(or neutral sputtering products, which would carry zero electric current) are of 
much less importance to overall electrostatic properties of the wake than the 
interplay between primary and secondary electron emission [Zimmerman et al. 
2011 ]. 



The crater depth her was chosen to be at least tens of Debye lengths Asw deep 
under the dynamic range of conditions imposed by the passing storm ( h cr and Aaw 
are measured in units of m). The condition Asw«h cr is representative of large 
craters at the lunar poles (A SM/ ~15 m for the typical solar wind); thus, a range of 
crater depths from 250-1000 m deep was simulated. The primary simulation data 
plotted in Figs. 2-5 and 7 were taken from the 250 m-deep crater subset; however, 
the effects of crater depth are discussed in Section 2.2 and shown in Fig. 6. The 
simulation grid spacing is chosen to be less than Asw during each CME stage. It is 
computationally unfeasible to vary the plasma conditions continuously over a 
simulated timescale of days, so each stage of the passing CME is simulated 
independently according to the solar wind parameters given in Table 1. In most 
cases, only a small amount of simulated time corresponding to tens of ion plasma 
periods is required to reach a quasi-steady state, after which the resultant fields and 
computed fluid quantities such as particle density and flow speed were time 
averaged for hundreds of ion plasma periods. 

2.2. Basic modulation of wake structure by passage of the storm 

Figure 2 shows basic results from the kinetic plasma simulations described 
in Section 2.1, including net charge concentration, electric field, and electric 
potential structure within a 250 m-deep crater during each stage of the passing 
storm. Under typical solar wind conditions (Fig 2a) electrons race into the crater 
ahead of the slower ions, forming a persistent non-neutral cloud of negative space 
charge. An ambipolar electric field sets up between the negative electron cloud and 
positive ions at the wake flank, and the thermal ion population enters the crater via 



acceleration by the ambipolar field. The entire ambipolar field region is roughly a 
few Asw deep. A negative surface potential forms all along the surface; however, the 
potential is most negative in the region where electrons reach the crater floor 
accompanied by relatively few protons. Farther downstream the plasma 
environment returns toward a "wall sheath" configuration in which the near-surface 
space charge density is positive, the surface potential is negative, and the horizontal 
variation in structure diminishes. It will be shown that the basic structural 
components of the wake including the ambipolar charge separation, ambipolar 
electric field and potential, and surface potential, as well as the flux of plasma 
particles to the surface, are modulated by changes in solar wind conditions as the 
storm passes by. 

It is worth noting that the system reaches a quasi-steady state when current 
balance (i.e. net zero electric current) is achieved everywhere along the surface. In 
the kinetic simulations employed, the particle currents at the surface are statistical 
in nature: at any given time and location an integer (or zero) number of incident 
particles perturbs the local charge distribution, emitting secondary electrons and 
generating a local plasma response. Thus, averaging over many snapshots of the 
particle flux to the surface produces a quasi-steady, statistical depiction of the 
surface fluxes, shown in Figure 3. Underneath the electron cloud the primary and 
secondary electron fluxes generally balance one another, being roughly an order of 
magnitude more intense than the local ion fluxes (except in the late cloud case 
where secondary emission is relatively low). Ions are also important in establishing 



current balance farther downstream because the plasma expansion delivers a 
significant population of ions as well as primary electrons to the surface. 

In the "shock" stage of the storm the bulk plasma density no is four times 
higher and the electron temperature T e is almost an order of magnitude higher than 
in the typical solar wind. The parameters no and T e have units of m -3 and K, 
respectively. It is easily seen from Fig. 2b that the resulting increase in electron flux 
into the crater strengthens the ambipolar field and significantly enhances the 
negative surface potential everywhere. Secondary emission, which is roughly linear 
in T e in the regimes considered herein, increases during the shock and enhances the 
charge density within the negative electron cloud region. In contrast with the other 
stages of the storm, there is a significant population of secondary electrons near the 
surface much farther downstream, as denoted in Fig 2b. In fact, the surface is 
charged positive because the average secondary yield is above unity for the 
enhanced primary electron energies of the present case. Interestingly, formation of 
the ambipolar potential causes the surface to float negatively with respect to the 
bulk plasma, although the surface continuously emits secondary electrons and 
remains positive. The emitted secondary population forms a cloud of negative 
charge near the surface but since the potential rises monotonically with increasing 
height all emitted secondary electrons escape the surface. The current due to 
secondary emission from the surface must be included when modeling charging of 
human exploration infrastructure (cf., Section 3). 

The "early cloud" stage follows the shock stage with plasma parameters 
similar to the typical solar wind. The resulting plasma wake is depicted in Fig. 2c, 



showing a dear resemblance to the typical solar wind stage with respect to the 
magnitude of the ambipolar charge separation and electric field and potential 
structure. Perhaps the most notable difference between the early cloud and typical 
solar wind cases is that the wake appears to be "stretched" horizontally during 
passage of the early cloud; this is due to an increase in solar wind speed relative 
to the characteristic ion vertical speed as represented by c s (measured in m S' 1 ). The 
issue of variations in spatial scale will be considered more comprehensively in 
Section 2.3. Finally, during the "late cloud" stage, no is relatively high (about a factor 
of 15 greater than the typical solar wind) and T e is relatively low (about a factor of 
four lower than the typical solar wind). The most notable effect is that Asw is 
significantly reduced, thereby compressing the vertical scale of the ambipolar region 
by nearly an order of magnitude. 

In all cases, the density profile just above the crater wall sets the initial 
conditions for plasma expansion into the crater. However, in the late cloud it is 
obvious that fine details of the plasma sheath on the upstream plateau (outside and 
upstream of the domains of Fig. 2) can propagate far downstream (e.g. the blue 
streamer of positive space charge in Fig. 2d). The same feature is evident to varying 
degrees in some of the other simulations, but these "positive streamers" have a 
negligible effect on the large-scale structure of their parent wakes. 

2.3. Dimensional analysis and detailed variation of wake structure 

There are numerous physical parameters that govern the structure of a 
plasma wake near a dielectric, charge-collecting surface. These parameters include 
those already defined, as well as the peak secondary electron emission coefficient So 



(dimensionless) and electron impact energy at peak emission Eo (eV), and are listed 
in Table 1. From this collection of physical variables four dimensionless ratios have 
been constructed: 
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where ks is Boltzmann’s constant (J kg' 1 K' 1 ) and e is the electron charge (C) which 
have been introduced to express the electron temperature in units of eV, 
commensurate with the units of Eo. The ratio a represents the characteristic inflow 
angle of the ions with respect to the solar wind drift velocity vector, and will be 
shown to govern the horizontal and vertical spatial scales of the wake (a is 
equivalently the ion Mach number). The ratio 7 captures the effective depth of the 
crater with respect to the thickness of the ambipolar charge separation region. The 
ratio £ characterizes the amount of secondary emission from the surface at a given 
solar wind temperature and constant secondary electron yield. It will be shown that 
together 7 and £ govern the electric field and potential structure of the wake by 
controlling the interplay between surface and ambipolar electric fields. The effects 
of do have been previously explored by Zimmerman et al., 2011 and are not further 
investigated herein. Values for the parameters (l)-(4) during the solar storm are 
given in Table 1. 



2.3.1 Variations in horizontal and vertical spatial scales of plasma wakes 

The overall spatial scale of the wake is governed by properties of the ions 
since they carry a large fraction of the total inertia in the system. As mentioned 
previously, the primary velocity scales governing ion inflow are the solar wind 
convection speed Vsw and the ion sound speed c s . The canonical theory of 
quasineutral plasma expansion [Farrell et al. 2010, Mora and Pellat 1979, and 
references therein] suggests that these velocity scales may be wrapped into a "self- 
similar" dimensionless ratio parameterizing the ion density in the expansion region, 

( 5 ) 

where x andy are the horizontal and vertical distance from the bottom corner of the 
crater wall (e.g. the top rim of the crater is located at x=0,y=h cr ). Note that v w /c s is 
simply the ion Mach number given by (1). For cold ions the number density may be 
written purely as a function of £ [Mora and Pellat, 1979], 



n i = n o ex p(£+l) 
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Expressions (5) and (6) are particularly useful because they capture the spatial 
scaling of the wake in terms of background solar wind flow conditions as well as the 
obstruction size. For instance, multiplying the velocities of (5) by the basic ion 
plasma timescale copr 1 -[vth e /A.De) 1/2 (units of s] yields the horizontal and vertical 
length scales of the system, which are respectively VswCOpf 1 and CsCOpf 1 , and are listed 
in Table 1 for the stages of the passing CME. In the present work CsCtyf 1 is equivalent 
to the plasma Debye length; however, this will not be the case in more general 
scenarios where the ions and electrons are not in thermal equilibrium. Varying 



either of the length scales v S w(Opf l or CsCOpf 1 (or, equivalently, the velocity scales Vsw or 
c s ) modifies the ratio (1) and stretches or compresses the geometry of the wake. 

In Figure 4 (top row) the normalization length scales Vswcopi' 1 and c s o>pf 1 are 
validated against the expanding ion density contours across the wide range of 
conditions encountered during passage of the CME. The non-dimensionalized 
density contours remain roughly unchanged during the entire storm, indicating that 
the chosen normalization length scales capture the overall spatial structure of the 
wake. Furthermore, the equal-density contours from theory (Eq. 6) are straight 
lines fanning out from the crater rim, and the simulated contours of ion density in 
Fig. 4 are indeed quite linear, at least downstream of x=v S w(Opf 1 . Upstream of this 
region the canonical theory cannot be applied because it does not account for the 
strong separation of charge near the crater wall; in this upstream region the ion 
density contours curve gently in the simulations. Finally, Fig. 4 (as well as the 
density contours of Fig. 2, top row) shows that within the expansion region the ion 
density falls off exponentially with decreasing height, which is consistent with 
expression (6). 

The mismatch between the density contours of the middle and top rows of 
Fig. 4 (and the top row of Fig. 2) fortifies the conclusion that electrons do not strictly 
follow the ion flow; the lighter, faster electrons not only outpace the ions entering 
the crater, but also produce and quickly respond to local electric fields. The low- 
energy bulk of the electron distribution penetrates a few /Lv into the wake and is 
reflected back into the solar wind by the ambipolar field, indicated by a "bulge” in 
the electron population below the wake flank. The fastest electrons escape the 



ambipolar region and continue toward the surface where they accumulate to form a 
surface electric field and associated negative potential. Fig. 4 (bottom row) 
demonstrates a topological similarity between the potential structure and the 
electron density, stemming from the fact that plasma electrons within the non- 
neutral region constitute the source term in Poisson's equation for electric potential. 
In short, while electrons play a critical role in driving the ambipolar expansion it is 
the properties of the ions - namely VswCtyf 1 and CsCOpf 1 - that determine the overall 
scale of the wake. Essentially, solar wind electron pressure generates the wake’s 
electric field against the structural backdrop provided by ion inertia. 

2.3.2 Variations in the electric potential and electric field of plasma wakes 

Figure 5a compares representative vertical profiles of the vertical electric 
field across the wide variety of conditions encountered during the CME. The 
profiles are taken at a horizontal location x=0 .SvswOfyf 1 to exploit the horizontal 
scaling of the wake, providing consistency across all cases (see Table 1 for values of 
VswtOpf 1 ). The peak ambipolar field Eamb is on the order of the characteristic electric 
field kBTe/eAsw, which is consistent with theory [Crow et al. 1975; Mora and Pellat 
1979]. At very low levels of secondary emission (as parameterized by £) the surface 
can accumulate a large net amount of negative charge from the primary electron 
cloud, and the resulting surface electric field E sr f grows to envelope the entire non- 
neutral region [cf. Zimmerman et al., 2011]. This effect is clearly evident in the case 
of the late-cloud, in which the nearly uniform E sr / extends from the surface into the 
ambipolar region. In other words, there are two components to the wake’s electric 



field: the ambipolar field and the surface field due to accumulated surface charge, 
and the parameter ^governs the balance between the two field components. 

In all of the cases herein, the ambipolar field is the strongest component, 
evidenced by Fig. 5a. Also, Fig. 5b shows the corresponding potential structure of 
each case, where it is clear that a total potential drop from plasma to surface on the 
order of tens of ksT e /e is possible. Within a few Asw below the wake flank the space 
potential 0 tends to reach a few keT e /e due to the ambipolar expansion process. 
During the late cloud stage the magnitude of 0 clearly continues to increase with 
depth into the crater. There is a combination of effects giving rise to the very large 
potential drop during the late cloud. First, T e falls considerably with respect to the 
other stages, causing secondary emission to drop via a reduction in £ and allowing 
E S rf to extend throughout the expansion region (e.g. the "flat" electric field profile 
deep within the crater). Second, the Debye length shrinks considerably, causing an 
increase in 77. As a result, Reintegrates over a larger depth with respect to the very 
shallow ambipolar region to form a large potential drop. 

For a larger collection of runs including craters of different depths, the effects 
of secondary emission and crater depth are captured by the two dimensionless 
ratios <^and 77. Fig. 6a shows that the normalized potential drop <p=e0tot/kBT e , or 
potential "amplification", responds to the amount of secondary emission as 
parameterized by £"[cf. Zimmerman et al. 2011]. In addition. Fig. 6b shows that q> 
generally becomes more negative with increasing crater depth. The simulations 
span two important regimes of 77 where the ambipolar region is either partially or 
fully formed. In the former case where the crater is effectively shallow (e.g., ?7<20) 



the vertical space required to form the full ambipolar charge separation (about 
lOAsw according to Crow et al. 1979) is hardly available. In this regime solar wind 
electrons must overcome E am b integrated over a few Asw; thus, the normalized 
potential <p is relatively small. In the opposite high -77 extreme (e.g. the late cloud), 
the crater is effectively much deeper than Asw, and even a weak E sr f can then 
integrate over many Asw to create a large <p. In the 77>100-200 regime the charging 
timescale within the electron cloud appears to be much longer than the simulated 
time, and the numerical limit on the electric potential (~307V) is approached 
asymptotically (cf. Zimmerman et al. 2011). Thus, the late cloud potentials shown in 
Fig. 6b are in a state of slowly approaching the simulation limit but it would take an 
unfeasible computation time for the potentials to saturate. However, when 
compared to the other stages of the CME, all of the late cloud potentials of Fig. 6b 
demonstrate that there is a clear increase in \<p\ with increasing crater depth as 
parameterized by 77 . 

3. Astronaut charging in the regional wake environment 

3.1. Object Charging Model 

The constituent particle fluxes from each stage, shown in Figure 3 and listed 
in detail in Table 2, were fed forward into an equivalent circuit model of astronaut 
charging to anticipate the electrical hazards of stepping across the crater surface. 
The charging model was originally developed in Jackson et al. 2011, and is modified 
herein to include secondary emission from the astronaut suit During each stage of 
the passing CME, the model provides an accounting of the triboelectric and plasma 
charges accumulated by an astronaut moving in periodic steps (every 10 seconds) 



about a given location on the crater floor. To find the object potential relative to the 
solar wind plasma potential, a charging equation is solved given the local particle 
fluxes and secondary emission coefficient, § e ff~ 0.5-0.99, for the astronaut suit The 
primary assumptions in the charging model are that the astronaut suit is smaller in 
scale than the local Debye length and is immersed in electron and ion currents 
consistent with the equilibrium surface charge state. The latter criterion is 
automatically satisfied by the particle fluxes everywhere along the surface in the PIC 
simulations, and a spacesuit of about 1.5 m in height is considered, which has been 
verified to be much smaller than the local Debye length throughout the plasma 
expansion region (cf. Table 1). 

3.2. Astronaut charging rates and existence of charge equilibrium 

Two spatial locations representing regional environmental extremes are 
chosen including an upstream location (xj=250 m) within the tenuous, non-neutral 
expansion region and a downstream location (x2=1500 m) where the plasma density 
and ion fluxes are significantly greater. The charging and discharging behavior of an 
astronaut stepping in the vicinity ofxj and X 2 is shown in Figure 7. In this set of runs, 
the astronaut suit is considered to be insulating, with a low effective "triboelectric 
capacitance”; that is, the boot sole voltage acquired via accumulation of triboelectric 
charge is -13.6 V/step, as in Jackson et al. 2011. The secondary emission coefficient 
of the insulating suit is 5 e g= 0.5 [Jackson et al. 2011, Farrell et al. 2010]. 

Under typical solar wind conditions within the expansion region (i.e., at xi) 
the suit potential becomes more negative with each step due to triboelectric charge 
accumulation (e.g. Fig. 7a). In the absence of strong plasma currents, triboelectric 



charge accumulation effectively becomes perpetual; that is, charge is not dissipated 
at a rate faster than the astronaut's 10 s-per-step cadence. Much farther 
downstream (i.e., at X 2 in the presence of enhanced particle fluxes including a 
significant contribution from solar wind ions) charge is more quickly dissipated and 
the astronaut will return to the same equilibrium potential of about -10 V, or 
roughly -T e , in the 10 s before taking the next step (e.g. Fig. 7b). 

Model results considering an astronaut exposed to the shock passage are 
shown in Figure 7c, d. Throughout the entire crater the astronaut's accumulated 
charge will quickly dissipate between individual steps due to the enhanced flux of all 
plasma species during this warm, dense stage. Overall, the astronaut's potential 
during the shock stage is more negative than in exposure to the typical solar wind 
because of the increased number of secondary electrons near the surface and the 
enhanced primary electron flux both due to increased no and T e . More specifically, 
the astronaut’s equilibrium potential is about -50 V if roving upstream atxj and -15 
V downstream at X 2 \ the difference is due primarily to an upstream surplus of 
primary and secondary plasma electrons relative to the local density of ions. 
Interestingly, Figure 3 shows that in the upstream electron cloud region the total 
electron flux surpasses the ion flux by roughly a factor of about 100 but farther 
downstream the ratio is only about 10; thus, when located upstream (underneath 
the negative electron cloud) the astronaut floats to a more negative potential. 

Figure 7e,f shows results for astronaut exposure to environmental conditions 
during the early cloud. The charging and dissipation characteristics are similar to 
those encountered in the presence of typical solar wind conditions, which is a 



consequence of the similarity in T e and no between the two stages. As in the typical 
solar wind case, the astronaut floats more negatively with each step when roving 
upstream and tribocharging becomes perpetual in the absence of adequate charge 
remediation by environmental currents. Farther downstream the astronaut reaches 
an equilibrium potential of about -10 V. 

Finally, Figure 7g,h shows that during the very tenuous, cold late cloud the 
astronaut can be expected to charge indefinitely while walking upstream in the 
vicinity of xi, similar to the typical solar wind and early cloud cases. When roving 
about location X 2 downstream, the astronaut is able to remediate charge and come 
to an equilibrium potential of about -10 V before taking the next step. 

The next set of runs considers the same insulating astronaut suit now 
exposed to a relatively high triboelectric capacitance scenario in which the boot sole 
potential changes by 136 V per step. The results are qualitatively similar to those of 
Figure 7. In general, at 10 steps per second the astronaut's triboelectric charge 
accumulation dominates over the environmental plasma currents, and the suit 
potential diverges, as in Fig. 7a. The average divergence rate is generally slowest in 
exposure to the downstream plasma, due to the increased plasma fluxes providing 
more significant charge remediation there. However, there is still an overall 
perpetual charging behavior irrespective of the plasma conditions. It is notable that 
in the present scenario of Fig. 7 the electric potential reaches significantly higher 
levels of -300 V to -3000 V during the same one-minute period simulated in Fig. 7. 
The only set of conditions under which the suit actually reaches charge equilibrium 
is within the downstream environment of the warm, dense shock where the particle 



fluxes are significantly enhanced; here, the equilibrium potential is about -50 V, as in 
the case of Fig. 7d where the solar wind current source dominates over the 
tribocharging process. In a real lunar wake astronaut charging would eventually be 
limited by some other factor such as electrical discharging via physical contact with 
or proximity to local infrastructure, the finite excursion time associated with a 
moonwalk, field emission from sharp edges on the suit, or nonlinear feedback 
limiting the local plasma fluxes. 

For comparison, the charging dynamics of a fully conducting suit were also 
considered, with a secondary emission coefficient <%= 0.99 characteristic of a 
conductor at the primary electron energies modeled herein [Jackson et al., 2011, 
Farrell et al., 2010]. The behavior closely resembles that of the insulating case 
therefore figures are not included. The principal difference between the conducting 
and insulating suits is in their composition relative to the surface [Jackson et al. 
2011], allowing greater triboelectric charge accumulation by the conducting suit 
counterbalanced by a simultaneous decrease in suit capacitance (cf. Jackson et al. 
2011 for more details). 

4. Implications for the volatile balance within a PSR 

The fact that solar wind protons are accelerated into the crater by the wake 
electric field has implications for the sequestration and loss of volatiles at the crater 
floor. It has previously been argued on the basis of laboratory experiments [e.g., 
Johnson, 1989] that for a crater floor nearly devoid of water ice the H 2 O sputtering 
yield induced by solar wind proton bombardment [at the 1 keV energies observed in 
the simulations) will be approximately 0.01 molecule per incident proton 



[Zimmerman et al. 2011, Farrell et al. 2010, Crider and Vondrak 2003, and 
references therein]; in this case, solar wind protons will represent a source of 
hydrogen to the surface. However, for a water-rich surface the sputtering yield for 
water increases to near unity [Johnson, 1989], so in this scenario proton 
bombardment would serve to scour water from the surface (e.g. Table 3). 

Other solar wind ions may play a role in sputtering a wide variety of volatiles 
during passage of the solar storm. Of particular interest in the present work is the 
increase in He ++ (doubly-ionized He) density to about 20% of the proton density 
during passage of the early cloud [Skoug et al. 1999]. Although He ++ ions were not 
explicitly included within the simulations, the self-similar formalism presented in 
Section 2.3 may be exploited to estimate the amount of sputtering induced by He ++ . 
The sputtering yield goes roughly as ( MZ 2 )* where the exponent s falls within the 
range 1-1.6 [Johnson, 1990; Elphic, 1991]; thus, the characteristic sputtering rate 
due to He ++ impacts is anticipated to be as high as (4-2 2 ) 1 6 ~85 times the rate due to 
proton impacts under typical solar wind conditions. By analogy with solar wind 
protons, the He++ impact rate will vary with position along the surface due to the 
ambipolar expansion process (e.g., Fig. 1). Assuming thermal equilibrium between 
plasma species, the He ++ sound speed (characterizing the He ++ inflow rate) will be 
half the proton sound speed; thus, the bulk population of He ++ ions will reach the 
surface at roughly twice the distance downstream. The characteristic He ++ flux 
during the early cloud is estimated to be about 0.1 times the characteristic proton 
flux - that is, (nHe++/no)(c S He++/c s )=(0.2)(0.5)=0.1 - giving an effective maximum 
sputtering yield as high as 85 (0.8 H20/proton)(0.1noc s )~7.3xl0 n H 2 O nr 2 s 1 for 



an icy regolith or 85 (0.01 H20/proton)(0.1noc s )~9.2xl0 9 H 2 O nr 2 s 1 for a pure 
silica surface. 

More generally, integrated sputtering over many extreme solar wind events 
could play a role in determining the abundance of hydrogen on equator-facing 
shadowed slopes (i.e., slopes with normal vector pointing toward the equatorial 
plane of the Moon); however, the average implantation and sputtering rates during 
a typical lunation are currently poorly constrained for real topography. 
Nevertheless, the present simulation results suggest that within narrow, deep 
topographic features any accumulated water and other volatile species will 
generally be shielded from sputtering by protons and heavier ions. A sub-crater 
within a larger PSR may provide the most favorable environment for volatile 
sequestration by combining the effects of topographic shielding with shielding from 
photons scattered at the rim of the parent crater, or "double thermal shielding” 
[Colaprete, 2012]. In this scenario the sub-crater floor is able to maintain an 
extremely low temperature at which accumulated water is maximally stable while 
also being protected from scouring by the highly-variable solar wind. Such a 
hypothesis could explain the relatively high abundance of water detected within the 
LCROSS impact plume [Colaprete et al., 2010] with respect to the water 
concentrations inferred from LEND epithermal neutron measurements [Mitrofanov 
et al., 2010]. Furthermore, a highly topography-sensitive process - such as regional 
plasma wake formation - influencing the sequestration of volatiles at the lunar poles 
could also help to explain the enigmatic lack of spatial correlation between the 
locations of highest inferred water concentrations and regions of permanent 



shadow [Mitrofanov et al., 2010], although a more detailed understanding of wake 
formation in the vicinity of complex lunar topography is needed. 

5. Conclusions 

It is found that the physical parameters affecting wake structure in the 
vicinity of a step-like lunar crater may be distilled into at least four leading 
dimensionless ratios including a=v S w/c s , r]=H/Asw, £=T e /Eo, and So, which represent 
the characteristic solar wind inflow angle, effective depth of the crater with respect 
to the ambipolar region, and secondary emission properties, respectively (Eq. 1-4). 
During a solar storm a stationary observer at the surface would observe changes in 
the particle fluxes and ambipolar field as the local self-similar parameter £ varies 
with changing v w , c Sl and A sw (Eq. 5-6). Solar wind electrons form and quickly 
respond to local electric fields within the structural backdrop provided by the 
expanding ion population, and secondary emission plays a prominent role in 
establishing equilibrium current balance in regions where electrons predominantly 
reach the surface. In particular, changing the constituent velocities v sw ,and c s within 
the ratio a affects the horizontal and vertical spatial scaling of the wake. Decreasing 
the ratio £ decreases the overall amount of secondary emission, and allows a very 
large negative potential to develop underneath the electron cloud, up to a 
magnitude of tens or hundreds of T e . Increasing the ratio 77 also tends to increase 
the wake potential drop. 

Environmental particle fluxes were fed forward from the kinetic simulations 
into a model of triboelectric astronaut charging. The astronaut is considered to be 
electrically isolated so that plasma and secondary emission currents provide the 



only paths to remediation of accumulated triboelectric charge. A stepping rate of 
one step per 10 s was assumed, with the astronaut gaining negative triboelectric 
charge at each step. An equilibrium suit potential exists only in cases with enhanced 
particle fluxes (through increased no or TV); if the plasma currents are relatively low 
then suit will collect triboelectric charge until some other process mitigates the 
significant electron buildup on the suit (e.g., field emission, secondary emission via 
cosmic rays, etc. [Zimmerman et al., 2011]). Thus, for an insulating suit the hot, 
dense shock stage of the CME is the most favorable for avoiding perpetual charging 
of the astronaut suit; however, this period is the most hazardous with respect to 
other factors such as high-energy particle fluxes and radiation. 

The wake formation process may play a leading role in volatile sequestration 
and sputtering within PSRs. Heavier, highly-charged ions such as He ++ observed 
within the early cloud stage of the CME by Skoug, 1999 would tend to scour water 
away from the surface. However, it is suggested that narrow, deep polar features or 
sub-craters within larger PSRs are prime environments for protecting accumulated 
volatiles via (1) topographic shielding from heavier ions and (2) double thermal 
shadowing in which scattered sunlight cannot reach the sub-crater surface to 
sublimate volatile species, Topographic modulation of solar wind particle fluxes 
may help to provide an explanation for the enigmatic spatial distribution of surficial 
water at the lunar poles; however, a more detailed understanding of wake formation 
downstream of real lunar topography, over many lunations and extreme solar wind 
events, is required to fully address this issue. Simulation efforts are underway to 



generalize the results presented herein to a broader range of topographic 
geometries and scales. 
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Figure Captions 

Figure 1 : Illustration of plasma expansion process including behavior of electrons, 

protons, and He ++ ions. 

Fi gure 2 : Simulated plasma conditions within a polar lunar crater during passage of 

the modeled CME. Plasma parameters during each stage of the storm are given in 



Table 1. Blue and red density contours correspond to logic proton and net electron 
density, respectively, each normalized to the bulk solar wind density [no) before 
applying the logarithm. *Net charge concentration levels in column (d) have been 
amplified by a factor of 10 to reveal the relatively weak charge separation along the 
wake flank. 

Fi gure 3 : Simulated particle fluxes at the surface of the crater of Figure 2. Dots 
represent the time-averaged flux at each surface point, while solid lines depict a 
running spatial average along the surface. Electrons and ions are moving toward the 
surface (positive absolute particle flux to the surface), while emitted secondary 
electrons are moving away from the surface (negative absolute particle flux). 

Fi gure 4 : Contours of logio ion density (top), logio electron density (middle), and 
electric potential (bottom) from the cases in Figure 2, with spatial units scaled by 
Oipi^Vsw and (O p f 1 c s along the horizontal and vertical axes, respectively. Particle 
densities are normalized by no before applying the logarithm, and electric potential 
is normalized to ksT e /e and plotted on a linear scale. 

Figure 5 : Vertical profiles of the normalized, vertical electric field and normalized 
electric potential taken atx=0.5<u p ,- 3 Vsw within the 2S0 m-deep crater of Figure 2. 

Fi gure 6 : (a) Normalized total electric potential drop at at x=0. 5 (Opi^Vsw versus 
amount of secondary electron emission for a collection of 250 m- ( 500 m-, and 1000 



m-deep craters. An increase in ^ 1 implies an effective decrease in secondary 
emission, (b) Normalized total electric potential drop atx=0.5w p rV<rw versus 
normalized crater depth. 


Fi gure 7 : Charging history of an insulating astronaut suit at xi=2 50 m (a,c,e,g) and 
*2=1500 m (b,d,f,h) within the 250 m-deep crater of Fig. 2, at an astronaut stepping 
rate of 1 step per 10 s. 


Tables 
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Early Cloud 

Late Cloud 
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20 
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Table 1 : Bulk plasma parameters, secondary emission parameters, horizontal and 
vertical scaling lengths, and leading dimensionless ratios during passage of the 
modeled solar storm. 
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Table 2 : Ion, electron, and secondary electron currents deep within the upstream 
wake expansion region (atxi=250m) and much farther downstream (atX2=1500m), 
taken from the kinetic simulations (cf. Fig. 3). 
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